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Abstract 

We consider a system of stochastic partial differential equations 
modeling heat conduction in a non-linear medium. We show global 
existence of solutions for the system in Sobolev spaces of low regularity, 
including spaces with norm beneath the energy norm. For the special 
case of thermal equilibrium, we also show the existence of an invariant 
measure (Gibbs state). 

1 Introduction 

In this article we consider the following system of partial differential equa- 
tions 

dt(j){x,t) = TT{x,t) 

dt7T{x,t) = {d'^ — l)(j){x,t) — fi(p'^{x,t) — r{t)a{x) 

dr{t) = {r{t)-{a,7r{t)))dt + V2fdu;{t). (1) 

In Eqs.(P) {(pjir) is a pair of scalar fields satisfying periodic boundary con- 
ditions with X £ [0, 27r]. The vector- valued functions a = (ai, • • • , ax) are 
given with aj G , for some 7 > 0, for all j. The vector r = (ri, • • • , r^) 
IS m R^, and CO = {uji, • • • , uk) is a ET-dimensional Brownian motion. The 
matrix T is given by T = diag(Ti, • • • , Tk), and Tj will be interpreted as a 
temperature. The parameter is a coupling constant; we will be primarily 
interested in the cases fi = (linear Klein-Gordon equation) and n > 
(nonlinear defocusing linear wave equation). 

The system of equations Q arises from a model for heat conduction in a 
nonlinear medium. It can be derived from first principles from a Hamiltonian 
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system which consists of K linear wave equations in R coupled to a nonlinear 
wave equation in [0,27r]. The total Hamiltonian is given by 



+ 



[0,27r] 




i + \cPix)\^ + |7r(x)|2) + i^ld^ixt dx 

+ E( / d::cUj{x)pj{x)dx] / dx4>{x)aj( 

~{ VJR / \J[0,27r] 

One assumes further that the initial conditions of the {uj,Vj), j = 1, ■ ■ ■ , K 
( "the reservoirs" ) are distributed according to Gibbs measures at tempera- 
tures Tj. These measures are (formally) expressed as 

Z"^exp(--^ / {\dxUj{x)\^ + \vj{x)\^) dxj ]]_ duj{x)dvj{x) , (3) 

and they are simply the product of a Wiener measure (for the fields uj) with 
a white noise measure (for the momenta Vj). 

We refer to [TJ or [201 for details on the derivation of equations (Q) 
from the Hamiltonian system @ with initial conditions (jSJ, at least in the 
case where the nonlinear wave equation is replaced by a chain of nonlinear 
oscillators (formally a discrete wave equation). In that case one obtains a 
set of stochastic ordinary differential equations. The derivation is essentially 
the same as for the model considered here and will not be repeated. We 
simply remark that the derivation of Markovian equations is possible due a 
particular choice of the coupling functions pj. 

In a series of papers [lIlIISlIiniEniEIlElinillHl about the chain of non- 
linear oscillators, the existence, uniqueness, and strong ergodic properties 
of invariant measures have been established. Moreover, a number of prop- 
erties of these invariant measures have been elucidated, such as existence of 
heat flow, positivity of entropy production, and symmetry properties of en- 
tropy production fluctuations. These invariant measures represent station- 
ary states which generalize Gibbs distributions to non-equilibrium situations 
where there is heat flow. Ultimately our goal is to establish similar proper- 
ties for the systems of equations . But we study here the more immediate 
problems of existence of global solutions- a prerequisite for studying the ex- 
istence of stationary states- and existence and invariance of an equilibrium 
(Gibbs) measure. 
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In the case of equilibrium, that is, when all temperatures are equal, 
Tj = T for all J = 1, • • • , K, we will prove below that there is an invariant 
state given formally by the (non-Gaussian) Gibbs measure 

dz. = Z-iexp f--L / ^\d,cP{x)\' + \cPix)\' + !^\c^ixr + \n{xr)dx] 
\ 2,1 7[0,27r] 2. j 

X exp {--^'T] dr J]^ d(f>{x)dTT{x) . (4) 

To make sense of this measure, one considers first the Gaussian measure 
for the case yu = 0. Its support is contained in H'^ x H^~^ x for 
any s < ^ and, with probability I, 4> is also a continuous function. Hence 
we can think of the measure ly as the measure which is absolutely contin- 
uous with respect to with a Radon-Nikodym derivative proportional to 
exp(— / |(/)(x)|^(ix/4r). We expect, but have by no means proved, that the 
invariant measure for different temperatures, if one exists, has similar sup- 
port properties. But with this intuition, it is appropriate to seek solutions 
of (P) in spaces of rough data x H'^^^ x R-^ with s < ^. Indeed we show 
the global existence of strong solutions, for 1/3 < s < 1 (see Corollarv 13.41 
and the remark following it). We believe that these spaces, with at least 
1/3 < s < 1/2, are natural to the invariant measure problem. 

Clearly, in these spaces no energy conservation (or bounds on the en- 
ergy growth/dissipation) is available. In recent years, however, Bourgain 
P, Keel and Tao jTH] and many others have developed techniques to show 
global existence for wave equations and other Hamiltonian PDE's in Sobolev 
spaces below the energy norm. A review of recent results with an extensive 
bibliography can be found in jB]. Here, we use and extend these methods 
to establish global existence of solutions for wave equations coupled to heat 
reservoirs, i.e. with noise and dissipation. 

In the last section, we show that solutions to an ultra-violet cut-off 
version of our system of equations, Eq.(^, converge as the cut-off is re- 
moved. This result is then applied to show that the equilibrium Gibbs state 
described above, Eq.ljl}, is indeed an invariant measure in the case of 
equilibrium. Note that Gibbs measures for nonlinear wave equations (and 
nonlinear Schrodinger equations) have been constructed and studied by sev- 
eral authors, (Lebowitz, Rose and Speer [TJ], Zhidkov (22|, McKean and 
Vaninsky 16^, Bourgain jSHH], Brydges and Slade [H]) but for isolated sys- 
tems only. Note that in these works Gibbs measures for any temperature 
are invariant while in our case the temperature is selected by the coupling 
to the reservoir. Our work is also related in spirit to various recent works 
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on the ergodic properties of randomly forced dissipative equations, see e.g. 
[H IZI 1101 1141 117j and others. The main and very important differences are 
that our equation is Hamiltonian rather than parabolic so that there is no 
intrinsic smoothing in the equations, and that the dissipation is very weak. 



1.1 Notation 

It is convenient to write our system as Bourgain does |lj. Set 

i 

u = (f> + —7r (5) 

where B is the operator defined B = \/—d'x + 1. Note that (/) = 5Rn and 
-Ivr = Qu are respectively the real and imaginary parts of u. Thus our 
differential equations can be written, 

1 „ 

idtU = Bu + —{^(p + ra) 
B 

dr{t) = -{r{t)-{a,-K{t)))dt + V2Tduj{t). (6) 

Let 

u(u;,t) = (M,r) = (</> + —7r,r), (7) 

and let Uo(a;,t) = {uo,ro){oj,t) be the corresponding solution to the differ- 
ential equations but with the non-linearity turned off, /U = 0. 
For a vector quantity u = (n, r), we introduce the norms 

Hull/,. = (r2 + ||n|||,.)'/'- (8) 
The energy of a vector u is defined by 

£{vi) = \\\ufH,+\r' + ^J{^ufdx. (9) 

2 Estimates for the Linear Wave Equation 

In this section we collect basic estimates for the linear system, /u = 0. These 
estimates actually establish global existence for this system. 

The first step is to consider the linear (dissipative) system without the 
random driving terms, 

duo . „ .1 

at B 

^ = {Ba,^Uo)-ro (10) 
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where SjUq is the imaginary part of Uq. Set Lq = Bq + P with 



-iB 
-1 





-i^\a) 




fill 



Symbohcally, the solution of this system Eas. H10|) is given by e*'^°Uo(0) with 
Lo = Bo + P. The system should be regarded as hnear in a function space 
of complex functions over the reals (so that 9 is linear). 

Lemma 2.1 Assume that a £ H'^ for some 7 > and < s < 1. For Aq 
sufficiently large depending on the a's only, (i?o + Ao)(io + Ao)*~"^(-Bo + Ao)~* 
acting m © R is defined as a bounded invertible operator. 

Proof: We have the following operator estimates (the operators acting in 



P- 



1 



P 



Bo + Xo + X 
1 



{Bo + Xo + X)iBo + Xoy 



< 



< 



c(Ao) 

(1 + A)7 

c(Ao) 
(1 + A)7' 



(12) 



with 7' = min{7 + s, 1} and with c(Ao) — > for Aq 00. By expanding the 
resolvent for Lq in a geometric series, convergent for c(Ao) < 1, one finds 
from these estimates that 



{Bo + Ao)(Lo + Ao + XyHBo + Xo)-' 
= {Bo + Xoy-'{Bo + Xo + X)-^ + 



c(Ao) 
(1 + A)- 



Using 



{Lo + Ao 



,s-l 



dX 



(13) 



(14) 



lo {Lo + Ao + A)Ai-^ 

with Cs a suitable normalizing constant, and integrating the previous equa- 
tion, we obtain 

{Bo + Xo){Lo + Xor-\Bo + Xo)-' = 1 + 0(c(Ao)), (15) 

which clearly is bounded. By choosing Ao large so that c(Ao) is sufficiently 
small, we see that {Bo + Xo){Lo + XoY-^{Bo + Aq)"^ is invertible. I 
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Lemma 2.2 Assume a £ H'^ , with 7>0, 0<s<l. There is a constant 
C3 depending only on s and the a 's, such that 

||e*^°u(0)||j^. <C3||u(0)||h» (16) 

for all time t. 
Proof: We have that 

£o{^^)^\(\\u\\l.+r^) (17) 

is a (degenerate) Liapunov function for the hnear system Eq. (|10|) . The 
lemma follows if we can show that for a suitably large constant Aq, £o{{Lo-\- 
Ao)*~^u) is equivalent to \\{Bo + Ao)*u||2, which is in turn equivalent to 
||u|||^i,. This is certainly the case if s = 1. For s < 1, this amounts to 
showing that [Bq + \q){Lo + Ao)*~^(-Bo + Ao)~** is a bounded invertible 
operator, which is the content of the previous lemma. ■ 

We now provide an estimate for the hnear stochastic evolution Uq. 

Lemma 2.3 Assume that 7>0, 0<,s<l and set ||uo(0)||//s = /3. There 
exist constants c and C , such that for A > c(5, 



(A - c/?)- 



P { sup ||u„(t)||^,. > A ^ < Cexp ( . (18) 



Remarks: The estimate is certainly not optimal. It does not account for the 
rapid dissipation of energy for small k modes of Uq. The lemma provides a 
global hound on the linear evolution, showing that it doesn't blow up, almost 
surely. 

Proof. The solution Uo(t) can be written in integral form, 

Uo{t) = re(*~*')^°Voda;(O + e*^°u„(0) 
Jo 

= woooit) + / e^*'''^^°LoWoU:{t')dt' + e*^°Uo(0) (19) 

JO 

by integration by parts, with Lg defined as in the deterministic case, Ea. Hllj) . 
and 

/ \ 
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By Lemma l2. 21 there is a constant C3 such that ||e*^°Uo(0)||//3 < C3/3, and 
||g(t-t )^ovo||//a < csllvollffs (which is finite) so that from the integral equa- 
tion Ea. H19() above, we obtain the estimate 



< \Wo\\HA^{t)\+C3t\\Lo^o\\H-SUp\uj{t')\+C3\\Uo{0)\\H^- (21) 

t'<t 

Thus we can write for a suitable constant c that 

\K{t)\\H^ <c{l + t)sup\uj{t')\+c(3, (22) 
t'<t 

with P the norm of the initial data Uo(0). 

Now if at some time t', with t' < t, we have that ||uo(t')||//s > A, then 

A-c/3 
c(l+i) 



evidently 3iffs < sup^/<j |w(t')l! ^-^d so, for A > c/3. 



p|sup||u„(0||h= > a| < pjsup|u;(Ol >4T-r7^1 
[t'<t J [t'<t c(l + t} J 

for yet another suitable constant C depending on the dimension of r. This 
concludes the proof of the lemma. ■ 

For later use, we also note here some simple Sobolev inequalities, all in 
one-dimension only. 

Lemma 2.4 For s > (1/2 — 1/p) and p > 2 there is a constant c = c{s,p) 
such that 

\\p<c\mH^. (24) 



Also, for < 9 < {1 — 1/p), there is a constant c such that 

||<Al|oo<c||<^||^(,„i)||<^||^-/. (25) 

Finally, for s > 1/6 and s' < min(0, 3s — 1), or s' = and s > 1/3, there is 
a constant c = c{s, s') such that for 0i, (f)2, 4>3 G , 

||010203|Ih^' < c||<?!)i||HH|02||HHI'/'3||i/''- (26) 
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Remark: The first inequality of the lemma actually holds with s = 1/3 and 
p = 6, as can be proved using the Hardy-Littlewood-Sobolev inequality. For 
convenience we will use this inequality as well, although it isn't essential 
for our purposes. But as a consequence of this remark, the last inequality 
Ineq.® holds for s' = 0, s = 1/3. 

Proof: The first inequality of the lemma is proved by estimating 

< ||(i + ^2)-.p'/2||^„||^||P;^ (27) 

with p' conjugate to p and r = 2/(2 — p'), all of which is bounded provided 
that sp'r > 1, i.e., s > i — ^. One then applies Hausdorff'- Young to obtain 
the first assertion of the lemma. 

The second inequality of the lemma is shown by first noting that 

^i-) = i E (i + J)i/2 (l + -')'/'^(-) (28) 
which by the Schwarz inequality gives the special case {9 = 0) 

uwoo <cmH^. (29) 

Also, we have that 

m\x)<p r\4>r'\cP'it)\dt + \4,\Piy). (30) 

Estimating the integral by ||i?5*|l2(^^i) ll</'ll_H'i ^^'^ then integrating this inequal- 
ity (|3fl|) with respect to y over [0, 27r], we get 

27r|0r(x) < 27rpmli;_,^mH^ + ll^llV-Dll-^lb, (31) 

so that 

||<^||oo<c||<^||J(-;_/f)||0||^/f. (32) 

The second inequality of the lemma is then obtained by interpolation be- 
tween Ineas. ()29l32|) . 

To prove the last inequality of the lemma, Ineq. (|26|) . we suppose each 
of the (pi^s is in with s > 1/6. Pick p' with s > ^ — ^ and, for later 

purposes, | < p' < |. By Ineq . (|?7|) above, each (f)i is in i^' , and the double 
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convolution 0i * 02 * </'3 is in i"^ for ^ = J? — 2 by Young's inequality. Note 

that 2 < r < cxd. It is then easy to check that (1 + n?Y l'^4>\ * (^2 * is in IP' 
provided s'r' < — 1 where p- + ^ = \. This is so if s' < (r' is positive) and 
s' < — = f — ^ = |r — |< 3(s + ^) — I = 3s — 1. The special case with 
s' = 0, s > 1/3 is an immediate consequence of the first inequality (j^lj). ■ 



3 Estimates for the Non-linear Equations 
3.1 Local Existence 

The Duhamel integral representation of the system equations for u, Eq.©, 
can be written 

w'j,,. ,33) 

Fix s with g < s < 1, and for > 1 let Vji{(3,t) be the set of functions 
defined 

VR{f3,t) ^ |u(-) GC([0,t],i7^)|||u(0)||H= </3 

and sup||u(t')||//= < -R/? ? , (34) 
t'<t J 

and let J-r{P, t) be the (probabilistic) event that the Duhamel integral equa- 
tion Ea. (|55|) has a unique strong solution in 'Dji{P, t). We have the following 
local existence result. 

Proposition 3.1 Assume ^ < s < 1. There exist constants c\, C2, C3 and 
C such that i/ u(0) satisfies \\u{0)\\h'' < P, R > Scs and t < ci/{R^0^), 
then 

P{W.t)} > 1 - Cexp • (35) 

Clearly, the sets J^ji{l3,t) are nested, J^R{P,t2) C .Fr(/3, ti) if ti < t2- 
The event J^r{P) = Un^R{P,t/n) is the event that u(-) exists for some pos- 
itive time, and, in this time, is no bigger than i?/?. An immediate corollary 
of the above proposition is that J-r{P) occurs with probability one. 

Corollary 3.2 For s > ^, local existence of the solution u(-) in holds 
almost surely, 

pi^Rm = 1. (36) 



10 



Proof of Proposition We have 



Jt-t')Lo 







< C3t/iSUp||03(t') 



t'<t 



1 



< c/xtsup||u(t')||H» < c^t{RI3y < -RP, 
t'<t <J 



(37) 



for a suitable constant c. Here, we have used Lemma 12.21 and the Sobolev 
inequality Ineq. (|26j) of Lemma l2.4l to estimate ||0^||j|^s-i, assuming that s > 
1/6 and using s — 1 < min(0,3s — 1). Also, we have chosen t < (3c/ii?^/3^)^^. 
Furthermore 





2Tduj{t') 



(38) 



using that the left side Inea . is bounded by c(l + t) supt/<( |<^(t')l fo'^ ^ 
suitable constant c (cf. the proof of Lemma l2.3|) : this condition holds with 
probability exceeding 



> 1-Cexp 



t{l + t) 



(39) 



for suitable constants C2, C. Finally, since R > 3c3 we have ||e*^°u(0)||H'' < 



C3/3 < Together with Ineqs. 



this implies that the right side of 



the Duhamel integral equation Ea. H33|) is a map of i) into itself. 

It remains to check that the right side of the Duhamel equation Eq. (|33)) 
is contractive for small t. But clearly for two functions ui, U2 G !?/?(/?, t), 

'2 respectively, and for ti < t, 

" -4m 





with real field parts 0i and 



dt' 







< 



csfit sup 
t'<t 



+ 0102 + 



02)(t')llH= 



< 3c/xt(i?/3)2sup||(ui -U2)(0IIh« 



(40) 



by Ineq. (|2(ij) of Lemma 12.41 for a suitable constant c. Thus the Duhamel 
integral is a contraction for t < l/(3c/xi2^/3^). 

In summary, if t < ci/ ^R?0^), for a suitably small constant ci, and if the 
stochastic integral estimate Inea. p8() holds, the right side of the Duhamel ex- 
pression maps I'ij(/3, t) into itself and it is a contraction, so that by the con- 
traction mapping theorem, Eq.® has a unique strong solution in I?r(/?, f). 
Ineq. (|38|l holds with probability at least that given in Ineq . (|39|) . | 
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3.2 Global Existence 

Finally, in this section we provide a global estimate for the non-linear stochas- 
tic evolution u{t). Following Bourgain's methods for the non-linear wave 
equation P, we set 

ujv(i) = {uN{t),f{t)) = {{u{t) - P>NUo{t)), (r - ro){t)), (41) 

with the positive integer N to be chosen later, P^n projection onto the 
Fourier modes {k : \k\ > N}. Here, u{t) and Uo{t) (the linear evolution) are 
assumed to begin with the same initial data, u(0) = Uo(0), and are driven 
by the same stochastic driving terms, so that they are not independent: 
they are coupled. The quantity /Ar(t) is defined as it would be in the pure 
deterministic case, 

lN{t) = £{uN{t)), (42) 
with £ the energy defined by equation @. Set 

0, = min{i(4s-l),^(l-s),7}. (43) 

We will assume below that a £ with 7 > 0, and that 1/3 < s < 1. 
Our main result is the following: 

Proposition 3.3 Let u(0) = Uo(0), with (5 = ||u(0)||h= = ||uo(0)||i^- . Fix 
> 1, > and 5 > so that 6 + 6 < 9^. There exist constants c, C and 
an No = No{R,l3) and r = t{R,0), such that if N > No and t < tN\ then 

p|sup/^(0 > ii/?2iv2(i-^)| < Cexp (44) 

Proof: The stochastic differential of I Nit) is given by 

dlN{t) = (^Q{BuN,IJ-i^uf - fi{mNf + roP<Na) -r'^ 

+ rQ{Buo,P<Na))dt. (45) 

In particular, there are no dm or du'^ = dt terms, hence the differential is 
the same as if we were just considering a wave equation with dissipation. 
We proceed to estimate the terms on the right side. 
We have that 

mufit) - {mNfm2 < c" {\\p>NUom2\\uN\L + \\p>NUoml) 

< C' {\\P>NUom2\\u\\TMHi + \\P>NUoit)\\l) 

< c [N-'\\uomHsiN{tf^^ + iKmh) (46) 



12 



for suitable constants d' ^c. Here we have used Sobolev inequalities, first 
Inea.()25|l. then Inea. 1)24(1 of Lemma 12.41 The other factors in Eg . (145(1 are 
readily estimated, and we get 

diNit) < c[N-'\\UomH^lN{tf^^ 

+{\ro{t)\+N'-^-'\\uomH^ + \\nomh)lN{t)'/') dt. (47) 



Now assuming that /Ar(t) < iJ/^^iV^Ci-^) and ||uo(t)||H= < N'^ , with 
< 9 < 6.^, = min{i(4s - 1), |(1 - s),^}, we obtain 

dlN{t) < Af2(i-s) (Are0(Ar|(i-4s)) ^ Ar3e0(^-(i-s)) ^ n^o{N-"')^ dt 

< N^(^-'^0{N^-'^')dt. (48) 

It follows that lN{t) < RfS^N^^^-"') for t < T with T = O(iV^), 6 kO^-O, 
provided that in this time interval, ||uo(t)||//s < . Said more precisely, 
given P, R, there exist an No{(3, R) and a r(J, 6, (3, R) such that for N > No 
lN{t) remains less than Rfi'^N'^^^^^^ for a time t, < t < tN^ , provided that 
||uo(t)||//s remains less than in this same time interval. 
Thus we have that 

p|sup/^(0 > i?/32iv2(i-^)| 

= P \ supliv(t') > Rfi'^N'^^^-'^ and sup \\Mo{t')\\H- > \ 
[t'<t t'<t J 

< p|sup||uo(t')||H= >iV'| 

for t < tN^ , by Inea. ((18(l of lemma 12.31 and appropriate new constants 
Ci and ci. After renaming of constants and taking No still larger so that 
No > 2c/3, the proof of the proposition is complete. ■ 

The proposition 13.31 and the lemma [2^ give us a global bound: 

Corollary 3.4 Let j3, 6 < 9^, R > 2 he fixed, as in the above Proposition. 

There exist constants, c, C and Ni = Ni((3,R,t), such that for any time t, 
and N > Ni, 

p|sup||u(t')||H^ > RPN'-^^ < Cexp (^--0^^ (50) 



13 



Remark: The same ideas, in particular estimating I Nit) for any N, can be 
used to prove global existence almost surely in the energy norm s = 1, but 
we do not write out precise statements here. 

Proof: We have that 

p\sup\\u{t')\\Hs> Rm^-'] < p\snp\\uN{t')\\Hi>lRf3N^-'\ 
[t'<t J [t'<t ^ J 

+p|sup||uo(t')||H« > Ji?/3iVi-4. (51) 
\t'<t ^ J 

Now the first probability on the right side is bounded by 

p{s.p/„(0 > lii=/J^iV^<.-)} < Cexp (-J^) (52) 

for N > No by Proposition 13.31 The second probability on the right side 
of Ineq. (|51|) is bounded by the estimate given in Lemma 12.31 with A = 
\R(5N^-'' » 0{N^). Thus this probability is negligible compared to the 
first term on the right side of Ineq. (|5H) . Enlarging C completes the proof of 
the corollary. ■ 



4 Large k Cut-ofF Systems and an Equilibrium In- 
variant Measure 

4.1 Convergence of Finite Dimensional Cut-off Systems 

We consider a cut-off version of our system Eqs.(^ where we retain Fourier 
modes {k} with \k\ < M, M a positive integer. Let UM{t) = {uM,fM){'t) = 
{(t>M + ■^T^MirM){t) denote a solution to the finite dimensional system 

idtUMit) = BuM{t) + ^P<M{fJ-(t>M{t) + rM{t)a) 
drM{t) = -{rM{t)-{P<Mar.^M{t)))dt + y/Tfduj{t) (53) 

for initial data um{^) £ P<A/-^^^-The solution vLuii) remains in P<^mL^ 
and is clearly in (for any s), since all Fourier coefficients UM,k = 
for I A; I > M. We remark that under the same assumptions on the coupling 
functions a, the conclusions of the previous section, Propositions 3.1. l3.3l and 
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Corollary 13.41 hold for solutions uniformly in M with respective initial 
data uj\,f(0) = P<Afu(0) for an initial u(0) € H'^, s > 1/3. In particular the 
arguments used there are equally valid for the cut-off systems. 

Fix s and let I'_r(/3, t) be the set of continuous functions defined in 
Eg ■ (|34|) ■ in particular functions {u} bounded in the i7*-norm by RP with 
||u(0)||//s < /?, and let GR{(3,t) be the probabilistic event defined 

GR{/3,t) = {u(-), um(-) G VR{p,t) for each M} (54) 

with u(-) the solution to Eq.(^ and um{-) the solution to Ea. ((^ . 

Proposition 4.1 Fix s > 1/3, a time t > 0, and So > s. Then {um{-)} 
converges strongly to u(-) inH^ uniformly on G ^{P , t)n{u \ ||u(0)||//so < /?}. 

Proof. For notational convenience we will replace the t of the Proposition 
statement by ti, and work on the time interval t £ [0,ti]. We will assume 
that in this time interval we have the a priori bounds ||u(t)||//s < RP and 
||um (i)||_f/« ^ RP, for all M and t. Given these bounds, we will actually show 
the stronger result that -P<mu — um — > strongly in and P>mu 
strongly in , uniformly in t and the initial data with ||u(0)||/fso < p. 
The quantities u and um satisfy the respective Duhamel relations 



and, for with initial data P<j\/u(0), 

We proceed to estimate (P<a/u — UM){t) in the //^-norm, for < f < ti. 

From the above integral formulae Eas. H55l56|) one sees that there will 
be contributions to (P<jv/u — UM){t) coming from integrals involving the 
non-linearity — (t)\^, (r — rM)ct, and (a, (tt — ttm))- We first consider the 
integral of the non-linearity and obtain a contribution to (P<a/u — UA/)(t) 



\{t-t')B. ( ^P<M{<t^^-(t>li)dt' 







(57) 

t 



<c{RPff \\P>Mu{t')\\Hs dt' + c{RPff \\{P<MU-UM){t')\\m dt' 

Jo Jo 
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for a suitable constant c by the Sobolev inequality Inea. ()26() with s' = 0, s > 
1/3. Next, the contribution to (P<mu— UA,/)(i) from the integral of (r— rA,/)a 
is bounded by 



,{t-t')Bo 



-g P<MCtdt' 





< c\\a\ 



\P<MVi-ViM){t')\\m dt'. 

(58) 

Finally, the contribution to (P<mu — ujv/)(i) from the integral of (a, (vr — 
tta/)) is bounded by 

ll^'>Ma||2 f' \\P>Mii{t')\\m dt' + f\~^'-''\a,BP:^MUo{t'))dt' 
Jo Jo 

+ \\a\\2 f \\{P<MU-VLM){t')\\mdt', (59) 

Jo 

where Uq is the solution to the linear equations (with dissipation and noise) , 
and u = u — Uq. The second of these integrals is bounded by 



|^'>Ma| 



Lo + l 



u(0) 



< c||P>Ma||2||u(0)||2 

< c\\P>Ma\\2P 



(60) 



by Lemma l2.ll fwith s = 0). Adding these contributions, Inea. H57l58l59() . 

and using the last estimate, we arrive at 

||(P<MU - UmMWh^ < c{a,RP) r ||(P<MU - UM){t')\\H^ dt' 

Jo 

+c{RPf [ \\P>Mu{t')\\Hs dt' + WP^Mah [ \\P>Muit')\\Hi dt' 

Jo Jo 



+cRp\\P>Mah, 



(61) 



where c{a,Rf3) is linear in ||a||2 and quadratic in i?/3. 

The inhomogeneous terms on the right side of this inequality (|6H1 . i.e., 
the second, third, and fourth terms, each go to zero, M — > oo uniformly 
in t and the data. Second term of i61\) : Consider the projection of the 
integral formula of Eq. (|55j) above onto P^mH'^. To control the non-linear 
contribution to ||P>j\/u(t')||//s, we use 



IP- 



\P 



>M9 llilfs"-! 



u 



(62) 
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by the Sobolev inequality H26|) . with s" chosen, s < s" < 1. The a- contri- 
bution presents httle difficulty and is 0{M~'^ R(5\\a\\H-i) , while the inhomo- 
geneous term is estimated ||P>Me*^°u(0)||Hs < C3M'*~''°||P>j\/u(0)||iyso < 
C3M^-^°/3, by Lemma lO Third term of We have that 



dt' 



P>M [-^nt ) + ^ " ) 

< f'\\P>M{u{t'f)\\2dt' + 2Rpt\\P>Mah. (63) 
Jo 



Here, P>M{u{t'f) = P>M((P<[M/3]u(t') + ^'>[M/3]u(i'))'), where [M/3] de- 
notes the greatest integer < M/3. Expanding this out, one sees that terms 
containing a factor Py[]\,f/^ju{t') go to zero uniformly for M ^ oo as in our 
analysis of the second term of l|HT|) . and the term P^m {P<[M/3]^{t')'^) is iden- 
tically zero. The Q-term is 0{M^''' Rf3\\a\\Hi)- It follows that ||P>Mu(t)||//i 
goes to uniformly. Fourth term of i61\) : This term, proportional to a, is 
also OiM-TRPWaWH^). 

Thus Inea. H61() . with each of its inhomogeneous terms going to zero uni- 
formly, M — > oo, implies via Gronwall's inequality that (P<mu — UM)(t) 
goes to zero in uniformly in t and the data. Since P>Mu(t) goes to zero 
uniformly in H'^ as we have seen above, we have that UMit) converges to u(t) 
uniformly in , provided that ||u(t)||//s and ||uj\//(i)ll-ff'' ^^ay less than RP 
for t < ti, i.e., are in GR{f3,ti), and the data u(0) satisfies ||u(0)||i^so < 



4.2 Equilibrium Invariant Measure 

We proceed now to show the existence of an invariant measure for the com- 
plete system Eqs.(^, but in equilibrium where all temperatures are equal 
to a common T. Let be the Gaussian measure referred to in the in- 
troduction, which at least formally is the invariant measure for the linear 
problem with the / (/)^-term turned off. With respect to v^, WuWjjs has finite 
expectation for s < 1/2, 

/ iMlsdu^ = 2r^(i + k'^y-^ < oo, (64) 

k 

hence ||<?!'||_H"'i is finite, By Sobolev inequality (|24() of Lemma 12.41 

U < c||(/)||hs for s > 1/4, so that as random variables ||</)||4 and ||P<m0||4 
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are also finite I/O- a.s. (Actually i^q is supported on continuous functions, 
but we do not need this here.) Set 



(65) 



with appropriate normalizations 

Zm = y exp (^-^ y" (Ai/fi^a:^ ^ = exp (^-^ ^ <f dx\j (66) 

and (pM = P<m4>- Since ||(/)||4 and ||P<Af</'||4 are finite a.s., the Radon- 
Nikodym factors exp (— ^ / (t>%i dx) and exp {—^ J (p^ dx) are bounded and 
positive a.s., and the normalizations Z and Zm are positive. The measures 
vm and v are absolutely continuous with respect to . 

The semigroup associated with the cut-off system Ea. (|53p acts invari- 
antly on functions {/(u)} of the form /(u) = g{u-M, ■ ■ ■ ,UM,1~), with g 
integrable. The measure I'm is an invariant measure for this semigroup, as 
can be checked by computing the generator of the process and showing that 
its adjoint annihilates i^m', we leave this exercise to the reader. 

We also have that limj\/^oo Zm = Z. This is the case by the bounded 
convergence theorem: Clearly the exponentials are bounded by one, and 

I WHi - \\4'Mh\ < {{(p - (pMh < c\\4> - (PmWh" ^ o, M ^ oo a.s. (67) 
so that 

exp (^-^ J 4>Mdx^ exp (^-^ J (j)'^ dx^ a.s. (68) 

Let / be a function in the norm closure X of functions depending con- 
tinuously on only a finite number of modes, 

X = Um{/ I / = g{u-M-, ■ ■ ■ , UM-if), g bounded continous}. (69) 

Then again we have by bounded convergence that 
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Thus, vm converges to in a weak-* sense. 

For later use, we also note a kind of tightness for the measures {i^m}; 
for s < 1/2, 

which is arbitrarily small for [3 large, uniformly in M, by Inea. H64|) . 
Finally, we address the invariance of v. Define the semigroups 

SV(u) ^ i^u[/(u(t))], 5i,/(u) ^ i?u[/(uM(t))], (73) 

where f € X. (We will assume here for definiteness that the Fourier modes 
^M,k{t) are simply constant in time for modes \k\ > M.) 

Proposition 4.2 (Equilibrium case) The measure v is invariant with re- 
spect to the semigroup in the sense that for f £ X , 

I S'fd,^ = I fd,^. (74) 

Proof: Choose s, 1/3 < s < 1/2, and let GR{(3,t) be the event defined in 
Eg . (1541) (using the H'^-norm). By Corollarv 13.41 we have that for an initial 
u with ||u||//s < P, 

P^iQ'kiP.t)} < Cexp ( (75) 

we identify the R here with 2A^^~* in the corollary statement (the R of the 
corollary being chosen equal to 2) and appropriately redefine the constant 
c. 

Now let So be chosen, with s < Sq < 1/2 and let /(u) be a bounded 
function continuous in the ff'^-norm of u. (Such functions are dense in X). 
For ||u||/fso < /3 and any e > 0, 

|S*/(u)-5i,/(u)| 



< 



E^[Xg,(p,t){f{<t)) - fiuMit)))] +2\\f\\MQR{P,t)] 



< e (76) 

for R and then M chosen sufficiently large, by the above probability estimate, 
and by the uniform convergence of Um to u on ^_r(/3, t)n{u | ||u(0)||h''o < /?}, 



19 



Proposition 14.11 Consequently, 5^/(u) 
in u, ||u||/fao < (3. 
Finally, 



5*/(u), for M — > oo uniformly 



S'fdu - J SiJduM 



< 



S'f{dv - dvM) 



+ / X{\\m\\hso<P} 



dvM 



dvM- 



|u!Ihso>/3} 



(77) 



The first term on the right side goes to zero by weak-* convergence of {z^a/} 
to I/, the last term can be made arbitrarily small for [5 suitably large by 
tightness Ineq. (|72|) . and the middle term then goes to zero by uniform con- 
vergence of S\,if for ||u||//so < /3. Thus 



S^fdv 



lim 



S\ifduM 



lim 

M^oo 



fdvM = / fdu 



(78) 



by the above Ineq . (f77|) . by invariance of um under (5*^/)*, and again by 
weak-* convergence of the {vm}- This completes the proof of invariance of 
u for functions depending continuously on u with respect to the //''-norm 
and, by density, invariance for all / € ■ 

Remarks: We emphasize that the question of ergodicity for this equilibrium 
measure v remains open, as does the existence of non-equilibrium invariant 
measures for differing temperatures. 
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